Maturational networks of human fetal brain activity reveal emerging connectivity patterns prior to ex-utero exposure.

A key feature of the fetal period is the rapid emergence of organised patterns of spontaneous brain activity. However, characterising this process in utero using functional MRI is inherently challenging and requires analytical methods which can capture the constituent developmental transformations. Here, we introduce a novel analytical framework, termed "maturational networks" (matnets), that achieves this by modelling functional networks as an emerging property of the developing brain. Compared to standard network analysis methods that assume consistent patterns of connectivity across development, our method incorporates age-related changes in connectivity directly into network estimation. We test its performance in a large neonatal sample, finding that the matnets approach characterises adult-like features of functional network architecture with a greater specificity than a standard group-ICA approach; for example, our approach is able to identify a nearly complete default mode network. In the in-utero brain, matnets enables us to reveal the richness of emerging functional connections and the hierarchy of their maturational relationships with remarkable anatomical specificity. We show that the associative areas play a central role within prenatal functional architecture, therefore indicating that functional connections of high-level associative areas start emerging prior to exposure to the extra-utero environment.

Structural and functional asymmetry of the neonatal cerebral cortex.

Features of brain asymmetry have been implicated in a broad range of cognitive processes; however, their origins are still poorly understood. Here we investigated cortical asymmetries in 442 healthy term-born neonates using structural and functional magnetic resonance images from the Developing Human Connectome Project. Our results demonstrate that the neonatal cortex is markedly asymmetric in both structure and function. Cortical asymmetries observed in the term cohort were contextualized in two ways: by comparing them against cortical asymmetries observed in 103 preterm neonates scanned at term-equivalent age, and by comparing structural asymmetries against those observed in 1,110 healthy young adults from the Human Connectome Project. While associations with preterm birth and biological sex were minimal, significant differences exist between birth and adulthood.

Amplitudes of resting-state functional networks - investigation into their correlates and biophysical properties.

Resting-state fMRI studies have shown that multiple functional networks, which consist of distributed brain regions that share synchronised spontaneous activity, co-exist in the brain. As these resting-state networks (RSNs) have been thought to reflect the brain's intrinsic functional organization, intersubject variability in the networks' spontaneous fluctuations may be associated with individuals' clinical, physiological, cognitive, and genetic traits. Here, we investigated resting-state fMRI data along with extensive clinical, lifestyle, and genetic data collected from 37,842 UK Biobank participants, with the object of elucidating intersubject variability in the fluctuation amplitudes of RSNs. Functional properties of the RSN amplitudes were first examined by analyzing correlations with the well-established between-network functional connectivity. It was found that a network amplitude is highly correlated with the mean strength of the functional connectivity that the network has with the other networks. Intersubject clustering analysis showed the amplitudes are most strongly correlated with age, cardiovascular factors, body composition, blood cell counts, lung function, and sex, with some differences in the correlation strengths between sensory and cognitive RSNs. Genome-wide association studies (GWASs) of RSN amplitudes identified several significant genetic variants reported in previous GWASs for their implications in sleep duration. We provide insight into key factors determining RSN amplitudes and demonstrate that intersubject variability of the amplitudes primarily originates from differences in temporal synchrony between functionally linked brain regions, rather than differences in the magnitude of raw voxelwise BOLD signal changes. This finding additionally revealed intriguing differences between sensory and cognitive RSNs with respect to sex effects on temporal synchrony and provided evidence suggesting that synchronous coactivations of functionally linked brain regions, and magnitudes of BOLD signal changes, may be related to different genetic mechanisms. These results underscore that intersubject variability of the amplitudes in health and disease need to be interpreted largely as a measure of the sum of within-network temporal synchrony and amplitudes of BOLD signals, with a dominant contribution from the former.

Development of neonatal brain functional centrality and alterations associated with preterm birth.

Formation of the functional connectome in early life underpins future learning and behavior. However, our understanding of how the functional organization of brain regions into interconnected hubs (centrality) matures in the early postnatal period is limited, especially in response to factors associated with adverse neurodevelopmental outcomes such as preterm birth. We characterized voxel-wise functional centrality (weighted degree) in 366 neonates from the Developing Human Connectome Project. We tested the hypothesis that functional centrality matures with age at scan in term-born babies and is disrupted by preterm birth. Finally, we asked whether neonatal functional centrality predicts general neurodevelopmental outcomes at 18 months. We report an age-related increase in functional centrality predominantly within visual regions and a decrease within the motor and auditory regions in term-born infants. Preterm-born infants scanned at term equivalent age had higher functional centrality predominantly within visual regions and lower measures in motor regions. Functional centrality was not related to outcome at 18 months old. Thus, preterm birth appears to affect functional centrality in regions undergoing substantial development during the perinatal period. Our work raises the question of whether these alterations are adaptive or disruptive and whether they predict neurodevelopmental characteristics that are more subtle or emerge later in life.

Integrating large-scale neuroimaging research datasets: Harmonisation of white matter hyperintensity measurements across Whitehall and UK Biobank datasets.

Large scale neuroimaging datasets present the possibility of providing normative distributions for a wide variety of neuroimaging markers, which would vastly improve the clinical utility of these measures. However, a major challenge is our current poor ability to integrate measures across different large-scale datasets, due to inconsistencies in imaging and non-imaging measures across the different protocols and populations. Here we explore the harmonisation of white matter hyperintensity (WMH) measures across two major studies of healthy elderly populations, the Whitehall II imaging sub-study and the UK Biobank. We identify pre-processing strategies that maximise the consistency across datasets and utilise multivariate regression to characterise study sample differences contributing to differences in WMH variations across studies. We also present a parser to harmonise WMH-relevant non-imaging variables across the two datasets. We show that we can provide highly calibrated WMH measures from these datasets with: (1) the inclusion of a number of specific standardised processing steps; and (2) appropriate modelling of sample differences through the alignment of demographic, cognitive and physiological variables. These results open up a wide range of applications for the study of WMHs and other neuroimaging markers across extensive databases of clinical data.

Quantifying noxious-evoked baseline sensitivity in neonates to optimise analgesic trials.

Despite the high burden of pain experienced by hospitalised neonates, there are few analgesics with proven efficacy. Testing analgesics in neonates is experimentally and ethically challenging and minimising the number of neonates required to demonstrate efficacy is essential. EEG (electroencephalography)-derived measures of noxious-evoked brain activity can be used to assess analgesic efficacy; however, as variability exists in neonate's responses to painful procedures, large sample sizes are often required. Here, we present an experimental paradigm to account for individual differences in noxious-evoked baseline sensitivity which can be used to improve the design of analgesic trials in neonates. The paradigm is developed and tested across four observational studies using clinical, experimental, and simulated data (92 neonates). We provide evidence of the efficacy of gentle brushing and paracetamol, substantiating the need for randomised controlled trials of these interventions. This work provides an important step towards safe, cost-effective clinical trials of analgesics in neonates.

Hospitalized newborns often undergo medical procedures, like blood tests, without pain relief. This can cause the baby to experience short-term distress that may have negative consequences later in life. However, testing the effects of pain relief in newborns is challenging because, unlike adults, they cannot report how much pain they are experiencing. One way to overcome this is to record the brain activity of newborns during a painful procedure and to see how these signals are modified following pain relief. Randomized controlled trials are the gold standard for these kinds of medical assessments, but require a high number of participants to account for individual differences in how babies respond to pain. Finding ways to reduce the size of pain control studies could lead to faster development of pain relief methods. Here, Cobo, Hartley et al. demonstrate a way to reduce the number of newborns needed to test potential pain-relieving interventions. In the experiments, the brain activity of nine babies was measured after a gentle poke and after a painful clinically required procedure. Cobo, Hartley et al. found that the babies' response to the gentle poke correlated with their response to pain. Further data analysis revealed that this information can be used to predict the variability in pain experienced by different newborns, reducing the number of participants needed for pain relief trials. Next, Cobo, Hartley et al. used this new approach in two pilot tests. One showed that gently stroking an infant's leg before blood is drawn from their heel reduced their brains' response to pain. The second showed that giving a baby the painkiller paracetamol lessened the brain's response to immunisation. The new approach identified by Cobo, Hartley et al. may enable smaller studies that can more quickly identify ways to reduce pain in babies. Furthermore, this work suggests that gentle brushing and paracetamol could provide pain relief for newborns undergoing hospital acute procedures. However, more formal clinical trials are needed to test the effectiveness of these two strategies.

The Developing Human Connectome Project: typical and disrupted perinatal functional connectivity.

The Developing Human Connectome Project is an Open Science project that provides the first large sample of neonatal functional MRI data with high temporal and spatial resolution. These data enable mapping of intrinsic functional connectivity between spatially distributed brain regions under normal and adverse perinatal circumstances, offering a framework to study the ontogeny of large-scale brain organization in humans. Here, we characterize in unprecedented detail the maturation and integrity of resting state networks (RSNs) at term-equivalent age in 337 infants (including 65 born preterm). First, we applied group independent component analysis to define 11 RSNs in term-born infants scanned at 43.5-44.5 weeks postmenstrual age (PMA). Adult-like topography was observed in RSNs encompassing primary sensorimotor, visual and auditory cortices. Among six higher-order, association RSNs, analogues of the adult networks for language and ocular control were identified, but a complete default mode network precursor was not. Next, we regressed the subject-level datasets from an independent cohort of infants scanned at 37-43.5 weeks PMA against the group-level RSNs to test for the effects of age, sex and preterm birth. Brain mapping in term-born infants revealed areas of positive association with age across four of six association RSNs, indicating active maturation in functional connectivity from 37 to 43.5 weeks PMA. Female infants showed increased connectivity in inferotemporal regions of the visual association network. Preterm birth was associated with striking impairments of functional connectivity across all RSNs in a dose-dependent manner; conversely, connectivity of the superior parietal lobules within the lateral motor network was abnormally increased in preterm infants, suggesting a possible mechanism for specific difficulties such as developmental coordination disorder, which occur frequently in preterm children. Overall, we found a robust, modular, symmetrical functional brain organization at normal term age. A complete set of adult-equivalent primary RSNs is already instated, alongside emerging connectivity in immature association RSNs, consistent with a primary-to-higher order ontogenetic sequence of brain development. The early developmental disruption imposed by preterm birth is associated with extensive alterations in functional connectivity.

Challenges and future directions for representations of functional brain organization.

A key principle of brain organization is the functional integration of brain regions into interconnected networks. Functional MRI scans acquired at rest offer insights into functional integration via patterns of coherent fluctuations in spontaneous activity, known as functional connectivity. These patterns have been studied intensively and have been linked to cognition and disease. However, the field is fractionated. Diverging analysis approaches have segregated the community into research silos, limiting the replication and clinical translation of findings. A primary source of this fractionation is the diversity of approaches used to reduce complex brain data into a lower-dimensional set of features for analysis and interpretation, which we refer to as brain representations. In this Primer, we provide an overview of different brain representations, lay out the challenges that have led to the fractionation of the field and that continue to form obstacles for convergence, and propose concrete guidelines to unite the field.

Inferring pain experience in infants using quantitative whole-brain functional MRI signatures: a cross-sectional, observational study.

Background: In the absence of verbal communication, it is challenging to infer an individual's sensory and emotional experience. In communicative adults, functional MRI (fMRI) has been used to develop multivariate brain activity signatures, which reliably capture elements of human pain experience. We aimed to translate whole-brain fMRI signatures that encode pain perception in adults to the newborn infant brain, to advance understanding of functional brain development and pain perception in early life. Methods: In this cross-sectional, observational study, we recruited adults at the University of Oxford (Oxford, UK) and infants on the postnatal wards of John Radcliffe Hospital (Oxford, UK). Healthy full-term infants were eligible for inclusion if they were clinically stable, self-ventilating in air, and had no neurological abnormalities. Infants were consecutively recruited in two cohorts (A and B) due to the installation of a new fMRI scanner using the same recruitment criteria. Adults (aged ≥18 years) were eligible if they were postgraduate students or staff at the University of Oxford. Participants were stimulated with low intensity nociceptive stimuli (64, 128, 256, and 512 mN in adults; 64 and 128 mN in infants) during acquisition of fMRI data. fMRI pain signatures (neurologic pain signature [NPS] and stimulus intensity independent pain signature-1 [SIIPS1]), and four control signatures (the vicarious pain signature, the picture-induced negative emotion signature [PINES], the social rejection signature, and a global signal signature) were applied directly to the adult data and translated to the infant brain. We assessed the concordance of the signatures with the brain responses of adults and infants using cosine similarity scores, and we assessed stimulus intensity encoding of the signature responses using a Spearman rank correlation test. We also assessed brain activity in pro-pain and anti-pain components of the signatures. Findings: Between May 22, 2013, and Jan 29, 2018, we recruited ten healthy participants to the adult cohort (five women and five men; mean age 28·3 years [range 23-36]), 15 infants to infant cohort A (six girls and nine boys; mean postnatal age 4 days [range 1-11]), and 22 infants to infant cohort B (11 girls and 11 boys; mean postnatal age 3 days [range 1-10]). The NPS was activated in both the adults and infants, and reliably encoded stimulus intensity. The NPS was activated in the adult cohort (p<0·0001) and both infant cohorts (p=0·048 for infant cohort A; p=0·001 for infant cohort B). The SIIPS1 was only expressed in adults. Pro-pain brain regions showed similar activation patterns in adults and infants, whereas responses in anti-pain brain regions were divergent. Interpretation: Basic intensity encoding of nociceptive information is similar in adults and infants. However, translation of adult brain signatures to infants indicated substantial differences in infant cerebral processing of nociceptive information, which might reflect their absence of expectation, motivation, and contextualisation associated with pain. This study expands the use of brain activity pain signatures to non-verbal patients and provides a potential research approach to assess the impact of analgesic interventions on brain function in infants. Funding: Wellcome Trust, Supporting the Sick Newborn and their Parents Medical Research Fund.

Modelling subject variability in the spatial and temporal characteristics of functional modes.

Recent work has highlighted the scale and ubiquity of subject variability in observations from functional MRI data (fMRI). Furthermore, it is highly likely that errors in the estimation of either the spatial presentation of, or the coupling between, functional regions can confound cross-subject analyses, making accurate and unbiased representations of functional data essential for interpreting any downstream analyses. Here, we extend the framework of probabilistic functional modes (PFMs) (Harrison et al., 2015) to capture cross-subject variability not only in the mode spatial maps, but also in the functional coupling between modes and in mode amplitudes. A new implementation of the inference now also allows for the analysis of modern, large-scale data sets, and the combined inference and analysis package, PROFUMO, is available from git.fmrib.ox.ac.uk/samh/profumo. A new implementation of the inference now also allows for the analysis of modern, large-scale data sets. Using simulated data, resting-state data from 1000 subjects collected as part of the Human Connectome Project (Van Essen et al., 2013), and an analysis of 14 subjects in a variety of continuous task-states (Kieliba et al., 2019), we demonstrate how PFMs are able to capture, within a single model, a rich description of how the spatio-temporal structure of resting-state fMRI activity varies across subjects. We also compare the new PFM model to the well established independent component analysis with dual regression (ICA-DR) pipeline. This reveals that, under PFM assumptions, much more of the (behaviorally relevant) cross-subject variability in fMRI activity should be attributed to the variability in spatial maps, and that, after accounting for this, functional coupling between modes primarily reflects current cognitive state. This has fundamental implications for the interpretation of cross-sectional studies of functional connectivity that do not capture cross-subject variability to the same extent as PFMs.

Large-scale intrinsic connectivity is consistent across varying task demands.

Measuring whole-brain functional connectivity patterns based on task-free ('resting-state') spontaneous fluctuations in the functional MRI (fMRI) signal is a standard approach to probing habitual brain states, independent of task-specific context. This view is supported by spatial correspondence between task- and rest-derived connectivity networks. Yet, it remains unclear whether intrinsic connectivity observed in a resting-state acquisition is persistent during task. Here, we sought to determine how changes in ongoing brain activation, elicited by task performance, impact the integrity of whole-brain functional connectivity patterns (commonly termed 'resting state networks'). We employed a 'steady-states' paradigm, in which participants continuously executed a specific task (without baseline periods). Participants underwent separate task-based (visual, motor and visuomotor) or task-free (resting) steady-state scans, each performed over a 5-minute period. This unique design allowed us to apply a set of traditional resting-state analyses to various task-states. In addition, a classical fMRI block-design was employed to identify individualized brain activation patterns for each task, allowing us to characterize how differing activation patterns across the steady-states impact whole-brain intrinsic connectivity patterns. By examining correlations across segregated brain regions (nodes) and the whole brain (using independent component analysis) using standard resting-state functional connectivity (FC) analysis, we show that the whole-brain network architecture characteristic of the resting-state is comparable across different steady-task states, despite striking inter-task changes in brain activation (signal amplitude). Changes in functional connectivity were detected locally, within the active networks. But to identify these local changes, the contributions of different FC networks to the global intrinsic connectivity pattern had to be isolated. Together, we show that intrinsic connectivity underlying the canonical resting-state networks is relatively stable even when participants are engaged in different tasks and is not limited to the resting-state.

Spatial parcellations, spectral filtering, and connectivity measures in fMRI: Optimizing for discrimination.

The analysis of Functional Connectivity (FC) is a key technique of fMRI, having been used to distinguish brain states and conditions. While many approaches to calculating FC are available, there have been few assessments of their differences, making it difficult to choose approaches, and compare results. Here, we assess the impact of methodological choices on discriminability, using a fully controlled data set of continuous active states involving basic visual and motor tasks, providing robust localized FC changes. We tested a range of anatomical and functional parcellations, including the AAL atlas, parcellations derived from the Human Connectome Project and Independent Component Analysis (ICA) of many dimensionalities. We measure amplitude, covariance, correlation, and regularized partial correlation under different temporal filtering choices. We evaluate features derived from these methods for discriminating states using MVPA. We find that multidimensional parcellations derived from functional data performed similarly, outperforming an anatomical atlas, with correlation and partial correlation (p < .05, FDR). Partial correlation, with appropriate regularization, outperformed correlation. Amplitude and covariance generally discriminated less well, although gave good results with high-dimensionality ICA. We found that discriminative FC properties are frequency specific; higher frequencies performed surprisingly well under certain configurations of atlas choices and dependency measures, with ICA-based parcellations revealing greater discriminability at high frequencies compared to other parcellations. Methodological choices in FC analyses can have a profound impact on results and can be selected to optimize accuracy, interpretability, and sharing of results. This work contributes to a basis for consistent selection of approaches to estimating and analyzing FC.

Artificial limb representation in amputees.

The human brain contains multiple hand-selective areas, in both the sensorimotor and visual systems. Could our brain repurpose neural resources, originally developed for supporting hand function, to represent and control artificial limbs? We studied individuals with congenital or acquired hand-loss (hereafter one-handers) using functional MRI. We show that the more one-handers use an artificial limb (prosthesis) in their everyday life, the stronger visual hand-selective areas in the lateral occipitotemporal cortex respond to prosthesis images. This was found even when one-handers were presented with images of active prostheses that share the functionality of the hand but not necessarily its visual features (e.g. a 'hook' prosthesis). Further, we show that daily prosthesis usage determines large-scale inter-network communication across hand-selective areas. This was demonstrated by increased resting state functional connectivity between visual and sensorimotor hand-selective areas, proportional to the intensiveness of everyday prosthesis usage. Further analysis revealed a 3-fold coupling between prosthesis activity, visuomotor connectivity and usage, suggesting a possible role for the motor system in shaping use-dependent representation in visual hand-selective areas, and/or vice versa. Moreover, able-bodied control participants who routinely observe prosthesis usage (albeit less intensively than the prosthesis users) showed significantly weaker associations between degree of prosthesis observation and visual cortex activity or connectivity. Together, our findings suggest that altered daily motor behaviour facilitates prosthesis-related visual processing and shapes communication across hand-selective areas. This neurophysiological substrate for prosthesis embodiment may inspire rehabilitation approaches to improve usage of existing substitutionary devices and aid implementation of future assistive and augmentative technologies.

Disambiguating brain functional connectivity.

Functional connectivity (FC) analyses of correlations of neural activity are used extensively in neuroimaging and electrophysiology to gain insights into neural interactions. However, analyses assessing changes in correlation fail to distinguish effects produced by sources as different as changes in neural signal amplitudes or noise levels. This ambiguity substantially diminishes the value of FC for inferring system properties and clinical states. Network modelling approaches may avoid ambiguities, but require specific assumptions. We present an enhancement to FC analysis with improved specificity of inferences, minimal assumptions and no reduction in flexibility. The Additive Signal Change (ASC) approach characterizes FC changes into certain prevalent classes of signal change that involve the input of additional signal to existing activity. With FMRI data, the approach reveals a rich diversity of signal changes underlying measured changes in FC, suggesting that it could clarify our current understanding of FC changes in many contexts. The ASC method can also be used to disambiguate other measures of dependency, such as regression and coherence, providing a flexible tool for the analysis of neural data.

The developing human connectome project: A minimal processing pipeline for neonatal cortical surface reconstruction.

The Developing Human Connectome Project (dHCP) seeks to create the first 4-dimensional connectome of early life. Understanding this connectome in detail may provide insights into normal as well as abnormal patterns of brain development. Following established best practices adopted by the WU-MINN Human Connectome Project (HCP), and pioneered by FreeSurfer, the project utilises cortical surface-based processing pipelines. In this paper, we propose a fully automated processing pipeline for the structural Magnetic Resonance Imaging (MRI) of the developing neonatal brain. This proposed pipeline consists of a refined framework for cortical and sub-cortical volume segmentation, cortical surface extraction, and cortical surface inflation, which has been specifically designed to address considerable differences between adult and neonatal brains, as imaged using MRI. Using the proposed pipeline our results demonstrate that images collected from 465 subjects ranging from 28 to 45 weeks post-menstrual age (PMA) can be processed fully automatically; generating cortical surface models that are topologically correct, and correspond well with manual evaluations of tissue boundaries in 85% of cases. Results improve on state-of-the-art neonatal tissue segmentation models and significant errors were found in only 2% of cases, where these corresponded to subjects with high motion. Downstream, these surfaces will enhance comparisons of functional and diffusion MRI datasets, supporting the modelling of emerging patterns of brain connectivity.

Low-threshold mechanoreceptors play a frequency-dependent dual role in subjective ratings of mechanical allodynia.

In the setting of injury, myelinated primary afferent fibers that normally signal light touch are thought to switch modality and instead signal pain. In the absence of injury, touch is perceived as more intense when firing rates of Aß afferents increase. However, it is not known if varying the firing rates of Aß afferents have any consequence to the perception of dynamic mechanical allodynia (DMA). We hypothesized that, in the setting of injury, the unpleasantness of DMA would be intensified as the firing rates of Aß afferents increase. Using a stimulus-response protocol established in normal skin, where an increase in brush velocity results in an increase of Aß afferent firing rates, we tested if brush velocity modulated the unpleasantness of capsaicin-induced DMA. We analyzed how changes in estimated low-threshold mechanoreceptor firing activity influenced perception and brain activity (functional MRI) of DMA. Brushing on normal skin was perceived as pleasant, but brushing on sensitized skin produced both painful and pleasant sensations. Surprisingly, there was an inverse relationship between Aß firing rates and unpleasantness such that brush stimuli that produced low firing rates were most painful and those that elicited high firing rates were rated as pleasant. Concurrently to this, we found increased cortical activity in response to low Aß firing rates in regions previously implicated in pain processing during brushing of sensitized skin, but not normal skin. We suggest that Aß signals do not merely switch modality to signal pain during injury. Instead, they exert a high- and low-frequency-dependent dual role in the injured state, with respectively both pleasant and unpleasant consequences. NEW & NOTEWORTHY We suggest that Aß signals do not simply switch modality to signal pain during injury but play a frequency-dependent and dual role in the injured state with both pleasant and unpleasant consequences. These results provide a framework to resolve the apparent paradox of how touch can inhibit pain, as proposed by the Gate Control Theory and the existence of dynamic mechanical allodynia.

Nociceptive brain activity as a measure of analgesic efficacy in infants.

Pain in infants is undertreated and poorly understood, representing a major clinical problem. In part, this is due to our inability to objectively measure pain in nonverbal populations. We present and validate an electroencephalography-based measure of infant nociceptive brain activity that is evoked by acute noxious stimulation and is sensitive to analgesic modulation. This measure should be valuable both for mechanistic investigations and for testing analgesic efficacy in the infant population.

Biomarkers, designs, and interpretations of resting-state fMRI in translational pharmacological research: A review of state-of-the-Art, challenges, and opportunities for studying brain chemistry.

A decade of research and development in resting-state functional MRI (RSfMRI) has opened new translational and clinical research frontiers. This review aims to bridge between technical and clinical researchers who seek reliable neuroimaging biomarkers for studying drug interactions with the brain. About 85 pharma-RSfMRI studies using BOLD signal (75% of all) or arterial spin labeling (ASL) were surveyed to investigate the acute effects of psychoactive drugs. Experimental designs and objectives include drug fingerprinting dose-response evaluation, biomarker validation and calibration, and translational studies. Common biomarkers in these studies include functional connectivity, graph metrics, cerebral blood flow and the amplitude and spectrum of BOLD fluctuations. Overall, RSfMRI-derived biomarkers seem to be sensitive to spatiotemporal dynamics of drug interactions with the brain. However, drugs cause both central and peripheral effects, thus exacerbate difficulties related to biological confounds, structured noise from motion and physiological confounds, as well as modeling and inference testing. Currently, these issues are not well explored, and heterogeneities in experimental design, data acquisition and preprocessing make comparative or meta-analysis of existing reports impossible. A unifying collaborative framework for data-sharing and data-mining is thus necessary for investigating the commonalities and differences in biomarker sensitivity and specificity, and establishing guidelines. Multimodal datasets including sham-placebo or active control sessions and repeated measurements of various psychometric, physiological, metabolic and neuroimaging phenotypes are essential for pharmacokinetic/pharmacodynamic modeling and interpretation of the findings. We provide a list of basic minimum and advanced options that can be considered in design and analyses of future pharma-RSfMRI studies. Hum Brain Mapp 38:2276-2325, 2017. © 2017 Wiley Periodicals, Inc.